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Abstract We investigate three-body recombination rates into deep dimers in cold atomic gases with 
large scattering length within hyper-spherical adiabatic zero-range approach. We derive closed analytic 
expressions for the rates for one- and two-species gases. Although the deep dimers are beyond the zero- 
range theory the latter can still describe the recombination into deep dimers by use of one additional 
short-range absorption parameter. The recombination rate, as function of the scattering length, retains 
the known universal behavior — the fourth power trend with characteristic log-periodic peaks 
however increasing the short-range absorption broadens the peaks until they are eventually completely 
smeared out. Increasing the heavy-to-light mass ratio in a two-species system decreases the distance 
between the peaks and increases the overal scale of the recombination rate. 


1 Introduction 

Three-body recombination is the principal source of loss of atoms in a trapped cold atomic gas. The 
loss rates are directly measurable as functions of the tunable interaction parameters. A number of 
recent experiments have provided results for various combinations of alkali atoms mm. 

The recombination reaction can conclude with either a shallow — weakly bound — dimer, or with 
a deep - strongly bound - dimer. The recombination rates into shallow dimers in cold atomic gases 
with large scattering length exhibit the well-known universal behavior - the rates, as function of the 
scattering length, show the fourth power trend with the characteristic log-periodic structures [51151 
111 ] . This behavior is universal - it can be described by a zero-range theory with only one parameter, 
the scattering length. The universal behavior is the consequence of the fact that for cold atoms with 
large scattering lengths both the three-body dynamics and the shallow dimers depend only on the 
long distance properties of the atomic interaction which is largely determined by the scattering length 
alone. 

The deep dimers, on the contrary, cannot be described by a zero-range — a long distance — theory 
as their properties depend on the short distance physics. For a long distance theory the recombination 
into deep dimers, which takes place at short distances, looks simply like a sink of probability at 
short distances. Therefore in zero-range theories the recombination into deep dimers can be described 
phenomenologically by introducing one extra parameter which creates shuch a sink of probability [51 

mm- 

We shall follow this philosophy — to introduce an extra phenomenological parameter to account 
for the loss of probability at short distances — and calculate the recombination rates as function of 
scattering length within a hyper-spherical adiabatic zero-range theory. 

What distinguishes us from previous investigations is that we derive a closed analytic expression 
which works for both one- and two-species cold gases. 
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2 Zero-range model 

2.1 Hyper-spherical adiabatic approximation 


Within the zero-range hyper-spherical adiabatic approximation the wave-function of the three- 
body system is written as 

&(p,Mp) = P~*f(p)$(p, Vp), (1) 

where the hyper-radius p and the hyper-angles Q p are defined in appendix[A] The hyper-radial function 
d>(p, flp) is normalized to one in hyper-angles and for large p loses its dependence on p and turns into 
a hyper-spherical harmonic. Its exact form is of no importance in the present context. 

The hyper-radial function f(p) satisfies the hyper-radial equation 


d 1 2 + v 2 (p) ~ \ 


2 inE 


dp- 


P A 


f{p) = 0 , 


( 2 ) 


where m is the mass-scale from the definition (311 of the hyper-radial coordinates, E is the energy 


of the three-body system, and v{p) is the angular eigenvalue given by the lowest root of the angular 
eigenvalue equation. The latter depends on the scattering lengths and masses of the particles. 

For a system of three identical particles — which corresponds to a one-species cold gas — the 
eigenvalue-equation is given as 
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where ai < 0 is the scattering leirgtlQ and p\ is dimensionless reduced mass from the definition (311 
of the hyper-spherical coordinate^ 

For a two-species cold gas we shall consider the common situation with one light particle with mass 
TO 3 and two identical heavy particles with masses mi = m 2 ■ We shall assume that only the scattering 
lengths between the light and the heavy particles, a\ = 02 , are large and that the scattering length 
a 3 between the heavy particles is, in comparison, negligibly small. In this case the angular eigenvalue 
equation is given as 


(vk\ 2 

— V COS —- . , -r 

V 2 / sm(2^i2) 


sin 


where 


4 > 12 = arctan 





( 5 ) 


The shape of the effective potential in the hyper-radial equation © for one- and two-species systems 
is illustrated on Fig. [l] it features an attractive pocket at short distances followed by a repulsive barrier 
at long distances - an archetypal shape, conducive to narrow resonances behind a barrier. 

The dividing point, po , between the two regions is located where the effective potential is equal 
zero, that is, given by the solution to the equation v{p$) = For the one-species system 


Po ~ 0.84|ai|. 


( 6 ) 


For two-species systems we find the following numerical approximation to po as function of the heavy- 
to-light mass ratio mi/m^ in the practically relevant region of mass ratios, 


Po 

M 


0.29 ( ^ 
m 3 


0.69 


( 7 ) 


The quality of the approximation is illustrated on Fig. [2j 


1 We use the atomic-physics sign convention where the absence of the shallow dimer corresponds to negative 
scattering length. 

2 For three identical particles we choose m = mi = m 2 = m 3 and, consequently, pi = A. 
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Fig. 2 The solution po to equation (v(po) 2 — \) = 0 for two-species system as function of the heavy-to-light 
mass ratio mi/mz together with a phenomenological fit. 


2.2 Asymptotic regions 


In the long-distance region, p> |ai|, the solution to eigenvalue-equations for both one- and two-species 
systems is v = 2. This leads to an asymptotic equation with a (half-integer) centrifugal barrier, 


dp 2 


4-1 


2 mE\ 

n 2 ) 


f(p) = o, 


( 8 ) 


with the solution 

P~^f{p) = M k P ) - tan (5)Y 2 (kp) , (9) 

where J 2 {np), ^(np) are the regular and irregular Bessel functions, Hk = \J2mE, and <5 is the phase- 
shift determined by the physics at short distances. 

In the short-distance region, p <C |ai|, the eigenvalue-equation ^ for the one-species system reduces 
to 

/ U7T \ 8 . (VK\ 

-^os( T ) + ^sm(-j=0, 
which has an imaginary root v = is where s ss 1.006. 


( 10 ) 
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Fig. 3 The square, i'(O ) 2 = —s 2 , of the solution to the short-distance eigenvalue equation (|ll| for two-species 
system as function of the heavy-to-light mass ratio mi/ m 3 , together with a phenomenological fit. 


The eigenvalue equation Q for two-species system in the short-distance region p <C |ai| reduces to 

(V 7T\ 2 

- v cos —- —7 -r 

V 2 / sm(2^i2) 

which also has an imaginary root v = is. However, the parameter s now depends on the heavy-to-light 
mass ratio 7711/7713 as illustrated on Fig. [3j increasing mass ratio increases s, and thus the attraction 
at short distances. 

Thus for both one- and two-species systems the hyper-radial equation at short distances acquires 
a super-attractive effective potential, 





f(p) = 0, 


( 12 ) 


where we have assumed that the energy is much smaller than the potential. The equation has two 
linearly independent solutions, 


1 J~PP ±; 


Vp e ±islnp , 


(13) 


or, equivalently, 


y/osin(slnp), v / pcos(s lnp) . 


(14) 


The three-body system with zero-range potentials must be regularized El at some short-range scale 
r. The regularized wave-function at p -C |oi| thus takes the form 


P 2 /( p ) 



(15) 


2.3 Qualitative solution to hyper-radial equation 


A qualitative solution to the hyper-radial equation can be obtained analytically by assuming that the 
two asymptotic solutions — being exact correspondingly at p> |ai| and p <C |cii| — are qualitatively 
correct also in the region p ~ |ai|. The phase-shift S can then be obtained by matching the two 
asymptotic solutions at po where attraction turns into repulsion. 

At p = po the short-distance function (15) arrives with the logarithmic derivative 


jfo h f) 

(P~®/) 


p—po 





(16) 
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where the factor /3 = Po/l a iI is given as 

f 0.84 


P = 


/ \ 0.69 

0.29 ( ^ j , for two-species systems . 


, for one-species system ; 


(17) 


This boundary condition causes the known log-periodic behavior of certain observables — as function 
of scattering length - in three-body systems with large scattering lengths. 

Matching the short-range and the long-range solutions at p = p 0 gives the equation to determine 
the phase-shift <5, 


J' 2 {np) ~ tan(c5 )Y^(np) 


cot s In 


/3\ai\ 


(18) 


J 2 (kp) — tan(5)Y 2 (Kp) /3|ai| \ r 

where the prime denotes the derivative of the Bessel function with respect to its argument. 

In the following we shall only be interested in the recombination rate at zero temperature, that is, 
at vanishing energies, k —> 0. In this regime the Bessel functions can be expanded as 


1 / 2 ( 2 ; —> 0) — —z 2 , Y 2 {z — > 0) — - z 2 . 

8 7 T 


and then the phase-shift in the low-energy limit is given as 

1 — /scot ^sln 


tan(<$) ss <5 ss (k/3|oi|) 


^ 1 + ^ s cot ^s In 


(19) 


( 20 ) 


3 Recombination rate 

3.1 Reaction rate in terms of hyper-spherical phase-shift 

The long-distance hyper-radial solution ([ 9 ]) in the limit p —> 00 becomes 

f{p -t 00) - 


2 g -iK.p+i(f> _|_ g e +inp-i4 


77 Kp 


where <j> is a real number and S = e 2lS . 

If S is a real number then IS) 2 = 1 and the hyper-radial flux density, 


3 =-1 


2 m 


di_dr 

dp dp 


f 


( 21 ) 


vanishes - there is no loss of probability if the phase-shift 6 is real. 

However, a complex phase shift, <5 = + i^sS, leads to a missing flux density, 


A i = 


hn 1 — |iS| 2 
m 27 TKp 


( 22 ) 


(23) 


where 

1- |S| 2 = 1 - e - 4SM « . (24) 

The total missing flux is obtained by integrating over the hyper-sphere. The integration element 
can be obtained from the relation 


(dVf 

dp 


1 


_ \ 3/2 d 3 Xid 3 yi 

H’iH'jk J 


dp 


E 5 


Ylrni 


zli 


p 5 df 2 j 


P ’ 


(25) 


where dV is the three-dimensional volume element and f df! p = 7r 3 . Besides this mass factors we also 
need to collect the (square of the) factor in front of J 2 (np) from the hyper-spherical expansion of the 
three-body plane-wave, normalized to unity within three-dimensional volume V, 


ik x x+ik„y _ 


V 


1 ( 2 k) 3 

V (up) 2 


J 2 (kp) + ... , 


(26) 















6 


where k 2 + k 2 = k 2 and x 2 + y 2 = p 2 . 

Collecting all factors gives the following expression for the probability loss A J per unit time — or 
recombination rate — within volume V, 


AJ = - 2 m* 


EmA 3/2 o5 2 h l- is 1 ! 2 

TT- ) ^ 7T--- 

[[rrii) m k 4 


(27) 


The factor V 2 AJ , often referred to as recombination constant K, can be written in terms of the 
three-body energy E = h 2 K 2 /(2m) as 


V 2 A J = m 3 


fYj m i 

\n*m 


3/2 


87r 2 


Ti 5 1- \S\ 2 

TO 3 E 2 


(28) 


directly determines the recombination coefficients in the density balance equations for cold gases. 


3.2 Recombination rate in zero-range model 


In the zero-range model, where deep dimers are absent, the description of recombination into deep 
dimers, unlike recombination into shallow dimer, requires one extra parameter to make the phase- 
shift S complex. This must be a short-range parameter as, indeed, to recombine into deep dimers 
all three atoms need to be at short-distances from each other - inside the range of their short-range 
potentials. This is also the region where the regularization of the zero-range three-body model takes 
place. Therefore the recombination range should be close to the regularization range r. The loss of 
probability — absorption — at short distances can be described by adding an imaginary part to 
the boundary condition at the regularization range. The easiest way to achieve this is to make the 
regularization parameter r complex, 

r —*■ re ~ ie . (29) 


The absorption parameter e determines the extent of the absorption. 

Inserting the complex regul ariz ation radius (291 into the expression for the phase-shift (20), t hen 
the phase-shift into expression (24) for 1— |5| 2 , and the latter, finally, into the recombination rate (28) 
gives the following expression for the recombination rate into deep dimers in our complex-boundary- 
condition zero-range model, 


V 2 AJ = m 3 


m i 

V TI m i 


3/2 


2 5 tt 2 —/3 4 |ai| 

in 


n 1 — cot In + ise) 

^ 1 + |scot ^sln + ise) 


(30) 


The expression has two phenomenological parameters - the regularization range r and the absorption 
intensity at short distances e. These parameters cannot be determined within the zero-range model 
and should be obtained by fitting the expression to experimental data. 

For one-species system the shapes of recombination rates as function of scattering length for different 
absorption parameters are shown on Fig. [4j For small absorption parameter the log-periodic peaks 
caused by the term sln(|ai|/r) factor are well formed. Increasing the absorption smears the peaks out 
until they eventually disappear completely. 

Two-species systems show similar smearing out of the log-periodic peaks with increasing the ab¬ 
sorption parameter. However, the shape of the rate is also influenced by the heavy-to-light mass ratio, 
as illustrated on Fig. [5] increasing the mass ratio decreases the distance between the peaks and also 
increases the overall scale of the rate. 
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Fig. 4 The recombination rate V 2 AJ |30| into deep dimers at zero temperature as function of scattering 
length a\ for different values of the absorption parameter t for one-species system. 



Fig. 5 The recombination rate V 2 AJ (|30| into deep dimers at zero temperature as function of scattering 
length asi for e = 0.1 and different values of the heavy-to-light mass ratio parameter m\/m,z for two-species 
systems. 


4 Conclusion 

We have investigated universal properties of the three-body recombination reaction rates into deep 
dimers within the zero-range hyper-spherical adiabatic approximation. Although a zero-range model 
cannot describe the deep dimers, it can still account for this reaction by introducing one extra param¬ 
eter, which describes short-range absorption — or loss of probability — in the hyper-radial equation. 

We have used a qualitative analytic solution of the hyper-radial equation comprised of the short- 
and long-range asymptotic solutions matched at the middle. This resulted in an analytic expression for 
the recombination rate as function of scattering length and the absorption parameter. We have given 
the expressions for both one- and two-species cold gases, the latter with one light and two identical 
heavy particles. 

For weak absorption the recombination rate as function of scattering length shows characteristic 
log-periodic peaks. Increasing the absorption broadens the peaks until they are eventually smeared 
out. 

For a meaningful comparison with experimental data the effects from the scattering length between 
the heavy particles need to be included [7J. The third scattering length modifies both the s and 















(3 parameters in the expression for the recombination rate in a non-trivial way through significant 
modification of the eigenvalue equaition. This should be a subject of a separate investigation. 


A Hyper-spherical coordinates 


Given the particle coordinates r; and masses mi, where i G {1,2,3}, the Jacobi coordinates are defined as 


Xi = - r fc ) , y; 


\/ l^jk I r i 


rrijVj + mt r k \ 

vrij +m k J 


where 

1 mjirik 1 + mt) 

— ; , fJ'jk — ; ; ■ 

m rrij + mt m m; + rrij + mt 

where m is an arbitrary mass-scale. 

The hyper-angular coordinates are the hyper-radius, p, defined as 

2 = 2 2 - r fc ) 1 2 3 4 5 6 7 8 9 10 11 12 

^ X ‘ m(mi + m 2 + m 3 ) 

and five hyper-angles, 

Q p = = arctan 

For identical particles the index can be omitted. 



(31) 


(32) 


(33) 


(34) 
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